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Abstract. Assuming the observed pulsar velocities to 
originate during asymmetric collapse of stellar cores, we 
compute the amplitude of gravitational waves emitted 
during type II and Ib supernova explosions and their de¬ 
tection rate from within a distance of 30 Mpc. At the 
rms-level of advanced laser interferometers h « 10“^^ at 
frequencies 300 — 1000 Hz the expected rate is about 1 per 
year. 
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Supernova explosions are among the most violent 
events known in nature. Supernovae type II and Ib are 
triggered by the gravitational collapse of evolved degen¬ 
erate cores of massive (> 10 Mq) stars and leave be¬ 
hind neutron stars or, possibly, black holes as a rem¬ 
nant. An enormous amount of energy is released during 
the collapse (roughly, the binding energy of a neutron 
star ^ 0.15Mqc^). Although most of this energy is car¬ 
ried away by neutrinos, a part of it can escape in the 
form of gravitational radiation provided that the collapse 
is not purely symmetrical. This collapse anisotropy (ex¬ 
pressed conveniently in terms of the fraction of the energy 
released, cMqc^) was first introduced by Ozernoy (1965) 
and Shklovskii (1970). 

It required more than 20 years for the idea of the 
collapse anisotropy to receive serious observational sup¬ 
port. First, the revision of the pulsar distance scale led 
to increasing the observed pulsar velocities by a factor of 
two up to values of about 400-500 km/s (Lyne & Lorimer 
1994). Observations of pulsars associated with young su¬ 
pernova remnants revealed even higher pulsar velocities 
up to 900 km/s (Frail et al. 1994). Such high velocities 
of young pulsars could hardly be obtained without addi¬ 
tional recoil (“kick velocity”) during the supernova explo¬ 
sion. Next, a direct evidence for a kick velocity of at least 
100 km/s has recently been obtained from observations of 


processing binary pulsar orbit in PSR J0045-7319 in the 
SMC (Kaspi et al. 1996). The collapse asymmetry result¬ 
ing in the kick velocity of a neutron star at birth may 
be due to different reasons. For example, recent calcula¬ 
tions of Burrows & Hayes (1996) has shown the ability 
of neutrino anisotropic emission to produce kick velocities 
of 400-500 km/s, as observed (see also Imshennik 1992; 
Bisnovatyi-Kogan 1993). 

One of the consequences of Lyne and Lorimer’s result 
was a recognition that the kick velocity imparted to a neu¬ 
tron star at birth has a power-law asymptotic form at high 
velocities (Lipunov, Postnov & Prokhorov 1996a,b). Using 
direct Monte-Carlo calculations of binary star evolution 
(the so-called “Scenario Machine”), they found that the 
Lyne-Lorimer pulsar transverse velocities are best repro¬ 
duced assuming the space (3-D) kick velocity distributed 
as 

2 , 0.19 

/(“'i “ (1+,,8.72)1/2 (1) 

where x = v/vq, Vq = 400 km/s. An important thing 
about having this distribution is that at high velocities 
(v > 500 km/s) it goes practically as that is much 

slower than the maxwellian tail oc exp(—u^), which some¬ 
times is assumed for the kick velocity distribution (e.g. 
Portegies Zwart & Spreeuw 1996). 

The anisotropic stellar collapse may be a source of 
gravitational waves (GW). To be detectable from the dis¬ 
tance of Virgo cluster ( 7 ^ 18 Mpc) by future ground- 
based laser interferometric and bar detectors, the GW 
energy emitted during a SN explosion should be about 
AE = 1O~^M0C^ (Thorne 1995; Schutz 1996). However, 
in numerical calculations this energy has been found to 
be very tiny, typically 10“® — 10“^ of the total energy 
output (see Muller 1996 for a review). Unfortunately, all 
these calculations are still far from being realistic con¬ 
sidering enormous numerical difficulties and uncertainties 
involved. 

Nevertheless, the recognition of high additional veloc¬ 
ities neutron stars acquire at birth suggests a way to esti- 
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mate e from observations, regardless of the unknown col¬ 
lapse anisotropy mechanism(s). Indeed, the kinetic energy 
of the neutron star motion, AE = where Mns 

is the neutron star mass, may be considered as a lower 
limit to the energy emitted in GW. Assuming Mj^s = 1-4 
Mq, we obtain e = 0.7 x {v/c)'^. With the distribution law 
(|^ the fraction of high velocity pulsars is 

poo 

P{>x)= (2) 

J X 

(the power-law asymptotics is valid for x > 2). For ex¬ 
ample, the fraction of pulsars with v > 1000 km/s is 
« 0.05. Of course, there should exist a maximum cut¬ 
off velocity; it must be higher than the maximum pulsar 
velocities observed, ^ 2000 km/s {x > 5), but its exact 
value only slightly changes the normalization coefficient. 
For e = 10“"^ {v ~ 3200 km/s) we obtain P(> 8) « 0.004, 
i.e. every 250-th neutron star may be born with the high 
velocity required. 

The neutron star galactic birth-rate is determined by 
that of massive (> 10 Mq) stars, which is ~ 1/25 years 
assuming Salpeter mass function f{M)dM oc 
and the mean galactic star formation rate of 1 Mq per 
year. Using the relations e = 0.7(u/c)^, x = i;/(400km/s) 
and Eq. (g) we may write galactic birth rate of high- 
velocity pulsars as 

R = 0.38/25 ~ 1-3 X lO-'^el^ '’® (3) 


(here and below subscripts indicate the quantities mea¬ 
sured in units of the corresponding power of ten, e.g. 
e_3 = e/10-3) 

The characteristic amplitude of the GW burst from a 
SN located at a distance r is (Thorne 1987; Schutz 1996) 


= 5 X 10 


-22 1/2 / IkHz^^^^ f 18Mpc^ 






(4) 


where fc is the characteristic frequency of the burst. This 
level should be detectable by the advanced LIGO inter¬ 
ferometer. The event rate from the volume V of the Uni¬ 
verse is (e.g. Phinney 1991; see also Lipunov, Postnov & 
Prokhorov 1996b) TZ 0.01 x (galactic rate) x (U/Mpc^). 
Using Eqs. (^) and (|^) we find 

= (lyr 

This means that we have chance to observe 1 GW-burst 
per year caused by SN II and Ib explosions at a level of 
hrms = 10“^^ provided that the energy is carried out at 
the frequency 300 Hz. 

Since he oc ~ v, the distribution of the GW-bursts 
caused by the SN explosions from a given distance r will 
have the form (j|). This enable us to compute the cumula¬ 
tive distribution of the number of events with the ampli¬ 
tude higher than a given one (log/V(> he) — log lie) using 



log ho 

Fig. 1. The logA(> he) — log/ic curve of the number of 
GW-bursts caused by core collapse supernovae as seen by a 
GW-detector with an rms-sensitivity he at the unitary sig- 
nal-to-noise ratio in 1-year integration time. The supernova 
rate are calculated for the baryonic matter distribution within 
the region of 30 Mpe from Sun according to Tully’s Nearby 
Galaxies Catalog (Tully 1988) with account for the SN-rate 
dependence on the morphological type of the galaxies. The 
distribution of SN-generated GW-amplitudes is taken assum¬ 
ing Lyne-Lorimer kick velocity distribution (see the text). The 
hatched region encompasses the assumed GW frequencies from 
1 kHz (lower boundary) to 300 Hz (upper boundary). The bro¬ 
ken curve represents the SN rate assuming constant GW energy 
fraction e = IQ—* (from Lipunov et al. 1995) with contributions 
of local group of galaxies indicated. 


the realistic baryon matter distribution within « 30 Mpe 
according to Tully’s Nearby Galaxies Gatalog (Tully 1988) 
(see Lipunov et al. 1995 for more detail). The result is pre¬ 
sented in Fig. 1. The hatched region corresponds to GW- 
frequencies of the signal from 300 Hz (upper boundary) to 
1 kHz (lower boundary). For comparison, we reproduce the 
log N — log h calculated for constant e = lO”"^ (Lipunov et 
al. 1995). Glearly, the more realistic e-distribution yields 
an order of magnitude smaller event rate because the rate 
of SN explosions with higher e strongly decreases. The 
shape of the curve is smoothed by the f{v) distribution. 
The mean e for Lyne-Lorimer velocity distribution (|^ is 
(e) oc / v'^f{v)dv « 4.4 x 10“®. 

It is seen from the figure that a few events per year are 
expected at the noise level of the advanced-LIGO rms- 
sensitivity ~ 10“^^. At the initial laser interferometers 
sensitivity level hrms ~ IQ-^^ expected SN rate is 
very low, ~ 0.01 per year. 
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In our analysis, we have not used any particular mech¬ 
anism for the collapse anisotropy and tried to rely on 
the observational properties of pulsars which are reliable 
products of supernova explosions. The only assumption 
we used is that the pulsar kick velocities are entirely due 
to the supernova explosion asymmetry. Of course, not ev¬ 
ery core collapse supernova may give rise to a pulsar. For 
example, some pulsar studies suggest the galactic pulsar 
birth-rate 1/125 — 1/250 yr“^ (Lorimer et al. 1993), much 
smaller than the galactic SN II rate. Additionally, some 
other mechanisms of the collapse anisotropy may oper¬ 
ate. This means that our results can be considered as a 
lower limit and the actual detection rate may be higher. 
We conclude that insofar as the observed pulsar velocities 
is a measure of the core collapse anisotropy, gravitational 
radiation bursts from supernova explosions should be de¬ 
tectable by the advanced laser interferometers and bar 
detectors in 1-year integration. 

The authors thank Prof. Vladimir Lipunov for helpful 
discussions. The work is supported by the INTAS grant 
No 93-3364, grant of Russian Fund for Basic Research No 
95-02-6053 and by the Center for Cosmoparticle Physics 
“COSMION” (Moscow, Russia). The work of S.Nazin was 
also supported by the grant of ISSEP No a96-1523. 

Note added in proof 

In a recent paper, Hartman (1996) argued that the ob¬ 
served transverse pulsar velocity distribution is repro¬ 
duced equally well assuming the Paczyhski (1990) distri¬ 
bution for initial pulsar velocity p{x)dx oc {l + x‘^)~‘^ with 
X = u/(600km/s). The same analysis as for our distri¬ 
bution (|^) shows that the resulting logA^(> he) — log he 
curve does not change appreciably giving appoximately 
the same (e) « 4 x 10“®. 
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